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ABSTRACT 

Gravitational lensing is generally treated in the geometric optics limit; however, when the 
wavelength of the radiation approaches or exceeds the Schwarzschild radius of the lens, 
diffraction becomes important. Although the magnification generated by diffractive grav- 
itational lensing is well understood, the astrometric signatures of diffractive microlensing 
are first derived in this paper along with a simple closed-form bound for the astrometric 
shift. This simple bound yields the maximal shifts for substellar lenses in solar neigh- 
bourhood observed at 20 GHz, accessible to high sensitivity, high angular resolution radio 
telescopes such as the proposed Square Kilometre Array (SKA). 
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1 INTRODUCTION 

Gravitational microlensing is a powerful tool to probe the con- 
stitu ents o f the solar neighbourhood, the Galaxy and b eyond 
fe.g. iWambsganss.2006, V In particular iGaudi fc Blooml (|2005i ') 
have propose astrometric microlensing as a technique t o de- 
tect su b-stellar objects in the solar neighbourhood, and iHeyl 
(|2010al lbl) argued that diffraction could provide important con- 
straints on lensing objects in the Kuiper belt and beyond. The 
combination of diffraction and astrometric lensing offers a new 
dimension to microlensing surveys. 

Several authors have examined g ravitational l ensing in- 
cluding the effects of diffractio n (e.g. Ing el fc Rubakha 1978; 



Elsteilll980l: iBontz fc Haugan|[l9 81: Dcgu chi fc Watsonl ll986!; 
Ulmer fc Goodmanll 19951 : iTakahashi, i2004i b However, the focus 



has almost entirely been on the magn if icatio n of the image. 
An exception is the work of iLabevrij l| 19941 ') that examines 
the possibility of using a planetary mass lens as a telescope. 
This letter will examine the astrometry of diffractive lensing; 
that is how does lensing affect the centroid of the light distri- 
bution including the effects of diffraction. As diffraction can 
amplify the magnification of a gravitational lens, so too does 
it increase the motion of the image. Measuring the motion of 
the image can provide constraints on the lens, source and their 
relative motion. 

The commissioning of the Square Kilometre Array (SKA) 
over the next decade will offer an unprecedented view of 



the radio sky. iKoopmans fc de Bruvnl (|2000l ) outlines some 
prospects for using the SKA to understand strongly lensed 
quasars and especially the small-scale structure of the lens- 
ing object. This letter also examines primarily the lensing of 
quasars but focuses on nearby lensing objects with the hopes 
to provide constraints on the number of small bodies in the 
solar neighbourhood. Such constraints are difficult to obtain 
otherwise. The letter is divided into a calculation (§ [2]) of the 
astrometric signature of lensing both in the diffractive and ge- 
ometric optics regimes, a description of the results (§ O and 
an evaluation of the prospects of observing this effect (§111). 



2 CALCULATIONS 



[Schneider etHI l|l992l ) give the magnification for a point 
source including diffraction 



l-if iu''/2 

u e ' 



Jo{uv)du 



(1) 



where u is the radial coordinate that integrates over the plane 
of the lens, Jo is the Bessel function of the first kind, and v is 
the impact parameter of the source relative to the lens. Both 
u and V are dimensionless and measure lengths in units of the 
reduced Fresnel length. 
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DdDd 



(2) 



where uid = (^(l + ^d) is the angular frequency of the radiation 
as it passes the lens. Hence the value of Ud which compares 
the angular size (rd/Dd) of the occulting portion of the lens 
to angular scale of its diffraction pattern {\/rd) is given by 



Ud = Td 



LOd 



c DdDds 



Ift 



(3) 



-Re 



(4) 



The parameter / is given by 

_ ojd Ds 
^ ^ c DdDds 

where Da is the distance to the source, Dd is the distance to 
the lens, Dds is the distance between the source and the lens, 
Zd is the redshift of the lens, and the Einstein radius is the 
characteristic length of the lens. 



Re = \ 2Rs 



DdDd 



Ds 



(5) 



for the Schwarzschild lens where Rs ~ 2GMd/(? where G is 
Newton's gravitational constant, Md is the mass of the lens 
and c is the speed of light; therefore, the value of 



(6) 



The parameter / compares the wavelength of the radiation 
to the Schwarzschild radius of the lens. The limit where the 
gravitational field of the lens is negligible is / = 0, so the effect 
of gravity on the form of the integral, Eq. U]), is quite modest. 
In a cosmological context all of the distances given are angular 
diameter distances. 

Here the occultation will be neglected (i.e. Ud — >■ 0). 
The integral can be calculated in closed form in terms of the 
confluent hypergeome t ric fu nction {\F\{a\b; z)) for Ud = 0. 
iGradshtevn fc Rvzhilj l| 19941 ) give relation (6.631.1) which in 
this particular case yields 

l-if T / X J 7r//4 jfir-/ In 2) /2ri /i ■f \ 

u 'e ' Jo(uv)du — e ' e ^ " i 2/ 

iF.(l-4;l;-4)- 
The result for / = is simply i exp(— 



2.1 Astrometry 

The gradient of the phase of the incoming radiation points to 
the apparent location of an unresolved source on the sky. This 
location on the image plane is given by 



ainJ 

— 

ov 



2~if iu-'/2 



Ji(uv)du 



(8) 



where denotes the imaginary part of the expression it pre- 
cedes. The first equality will also hold for an asymmetric lens 



where (u) 
holds 



-SiVlny. If = the following expression 



iFi 2 



V 2) ,Fr(l-il-l--i^ 



(9) 



where for f — the quantity in the brackets is unity and 5R 
denotes the real part of the expression in brackets. 

For values of /, w ^ 1 the ratio of the hypergeometric 
functions can b e conveniently ap proximated by Gauss's con- 
tinued fraction (|Cuvt et al.ll2008l ') 



iFi{a + l;b+l- 
iFi (a; b; z) 



.(10) 



b{b+l) 



a + 1 



1 + 



(fe + l)(6 + 2) 



b-a + 1 
{b + 2){b + ?,f 

a + 2 
, (& + 3)(6 + 4)' 



Although several techniques exist to deter mine the range 
of th e confluent hypergeometric function (e.g. iKarp fc SitnikI 
|2009| ). it is simpler in this case to resort to numerical experi- 
mentation, it appears that the follow inequality obtains 



< u (u - i;) s; 2/ 



(11) 



with the value oscillating between the two extremes. This 
yields useful estimates for the magnitude of the astrometric 
shift from diffractive lensing. Furthermore, the value in the 
brackets of Eq. [5] is purely real at these extrema, so they are 
also extrema of the magnification. 

2.2 Physical Optics 

For values of /, u ^ 1 the value given by Eq. ([8]) may be 
estimated by using a physical optics approximation. In partic- 
ular the square roo t of Eq. p]) may be a pproximated up to a 
constant phase by l|Schneider et al.lll992l ) 



I ^ e^^^lJ.^ 
where 



+ e 



t/2) 



U± 



and 



i (i; ± ^v^ + 4/ 
/ In \ u± 



ul 



u±v 



u± 

V 



dv 



+ 2/ 



± 1 



(12) 

(13) 
(14) 

(15) 



The values u± are defined to be positive and negative while 
IJ,± is always positive. The two images have opposite parity. 
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It is natural to interpret this as a negative value of the mag- 
nification for one of the images; hence there is an additional 
7r/2 of geometric phase of the positive term relative to the 
negative term. This choice may seem rather arbitrary, but it 
results from the stationary phase approximation of Eq. ^ 
where one extremum (m-) is a saddle point. This allows a 
simple expression like Eq. (I12p to approximate the results of 
Eq. IT]) accurately for large values of /. For more complicated 
lens geometries, th e phase lag is propor tional to the Morse 
index of the image (|Schneider et al.lll992l l. 

Furthermore, essentially by design the following holds 
d(f)±/dv = ~u±. Combining these results with Eq. ((S]) yields 
an estimate for 



+ fl-U- + 



(16) 



\l/2 



(u+ + U-) cos 



1 / ci In /i+ d In /!+ 



dv 



dv 



sin A(j!) 



where at this level of approximation the total magnification is 
+ 2/(1+ cos A(?!)) 



and 



+ — and fi+fi- 



„2 (j,2 ^ 4J) ■ 



(17) 



(18) 



The various definitions allow some further simplifications 
yielding 



■u +f{S + cos A0) ^ 2/ sin A(j> 



om 

4 y 

V jv'' + 3/) 
i;2 + 2/ 



V (1)2 + 4/) 



(19) 



The result from geometric optics obtains by neglecting the 
terms with A(j) yielding, 



(20) 



The maximum displacement due to lensing in the geometric 
limit is S-^/^i^is at « = 2^^^Re. 



3 RESULTS 

Diffractive effects can have a dramatic effect on the trajectory 
of images of gravitationally lensed sources. In particular from 
Fig [T] is it apparent that the maximal displacement is much 
larger when diffractive effects are considered. As in the geo- 
metric limit, the centroid lies along the line connecting the 
centre of the lens and the source. Furthermore, the centroid 
lies further from the centre of the lens than the source. The ob- 
served oscillations point back to the location of the lens, so the 
detection of three oscillations combined with the presumably 




Figure 1. The figure depicts the paths of images and sources for 
various impact parameters with the angular position of the lens 
fixed. The dashed circle denotes the Einstein radius of the lens. 
The red line gives the path of the source (unlensed). The green 
curve gives the path of the image without diffractive effects. The 
black curve gives the path of the image with diffraction for (/ = 10). 
The impact parameters from top to bottom are 3Re, Ret0.5Re 
and 4Re ■ The dashed red line indicates how to translate a location 
along the path of the source to the location of the image centroid. 



known proper motion of the source determines the impact pa- 
rameter between the source and lens, the proper motion, mass 
and distance of the lens unequivocally. 

Fig. El shows that the displacement of the image centroid 
from the source location oscillates between no displacement 
and 2R%/v outward. Furthermore, the minimal displacement 
occurs at a maxima in the magnification. The maxima of 
V {u — v) (where the black curves touch the blue curve) oc- 
cur at a minima of the magnification. In particular because for 
small values of v and large values of / the magnification is w ell 
approximated by a Bessel function l|Schneider et al.lll992l ). it 
is straightforward to estimate the peak displacement that oc- 
curs near the first zero of the Bessel function Jo (a;) at x m 2.4 
to be 



{n- 



0.83fRE at 



Re 



2.4/"^ for / > 1. 



(21) 



For smaller values of /, the peak displacement occurs for 
smaller values of v than given by this formula, and therefore 
the displacement is larger than given here. In particular, the 
displacement is larger than the geometric limit (Eq. I20p for 
/ > 0.17. The maxima for / — 0.17 occurs at v/Re ~ 5.5 as 
opposed to 14 as estimated from Eq. (|2ip . 
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Figure 2. The total magnification (red) and the difference (black) 
between the apparent radial position of the image centroid and the 
location of the source as a function of source location relative to the 
centre of the lens. The green curve gives the displacement according 
to geometric optics. The rapid oscillations are for / = 10 as in Fig.[l] 
and the slower oscillations are for / = 1. The blue curve gives the 
envelope of 2R^/v. 



The envelope of the displacement, 2R%/v, is robust re- 
gardless of the value of / or v, so it is natural to focus on the 
displacement by dividing by the size of the envelope, yield- 
ing Fig. [3l The approximation from physical optics is depicted 
by the dashed curve and follows the accurate calculation for 
/ = 10 very closely. However, for / = 1 the agreement is 
much poorer. Furthermore, the physical optics approximation 
does not precisely follow the simple envelope as the deviations 
below zero and above two manifest. Even if one is more care- 
ful and approximates the magnification as a Bessel function, 
the envelope only obtains approximately. The presence of the 
strict envelope results from an apparently thus-far unknown 
property of the hypergeometric functions and allows useful 
approximation of the possible signal. 

Fig. [3] also shows the effect of a finite source size to wash 
out the observed oscillations toward the geometric optics re- 
sult. It is not surprising that the rapid oscillations suffer a 
greater decrement than for / = 1. However, observational real- 
ities push the use of higher frequency observations to get finer 
angular resolution. For a fixed source size and impact param- 
eter, the size of the oscillations is proportional to f~^^^ while 
the angular resolution of a given telescope is proportional to 
frequency, increasing with /. Consequently if angular resolu- 
tion is the only factor in the accuracy of determining the cen- 




12 3 4 

v/Re 



Figure 3. The difference between the apparent radial position of 
the image centroid and the location of the source as a function of 
source location relative to the centre of the lens. The result for 
/ = 1 is the slowly varying sinusoidal curve and / = 10 is the 
more rapidly varying one. The dashed curves give the results using 
the physical optics approximation from Eq. I I19I I. For / = 10 it 
is nearly indistinguishable from complete result. The displacement 
from geometric optics is plotted in green. Notice for / = 10 there 
are about three peaks over a length of one Einstein radius for small 
values of v and more for larger values. The other colours assume 
that the angular radius of the source equals the Einstein radius. 
Blue is the geometric optics result, red is for / = 10 and magenta 
is / = 1. 



troid, the signal-to-noise of the measurement of the centroid 
is proportional to f~^^^; it makes more sense to perform the 
measurement at lower frequency. On the other hand with in- 
creased flux, the centroid can be determined more accurately, 
so this conclusion could change depending on the spectrum of 
the object. 

Determining the centroid of an object's emission is gen- 
erally more difficult than measuring the flux itself; therefore, 
searching for the diffractive flux variation would generally be 
more fruitful than looking for an astrometric signature, unless 
the flux from the source is inherently noisy making the os - 
cillation in the magniflcation difficult to detect. iHevll l|2010bl ) 
outlines using quasars as powerful tools to detect diffractive 
microlensing. Flat spectrum radio quasars generally have high 
brightness temperatures, so large fluxes from small solid an- 
gles. This can dramatically increase the expected signal-to- 
noise ratio for a diffractive microlensing event. On the other 
hand the flux from the quasar may b e inherently n oisy domi- 
nating the detector noise upon which iHevll (|2010bl ) focus. For 
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such objects the astrometric signature of diffractive microlens- 
ing is a powerful tool. 

Because of the envelope of the oscillation, it is straightfor- 
ward to estimate the magnitude of the displacement oscillation 
given the impact parameter of the source relative to the lens 
(6) and the properties of the lens itself. In particular for a 
planetary mass lens the following obtains 



{u — v) 




as along as the source may be considered compact compared 
to the diffraction fringes. In particular the SKA is expected to 
have an angular resoluti on of about 5 mas at 20 GHz, or / ~ 7 
for an Earth-mass lens (jSchilizzi et al.ll2007l '). Whether or not 
SKA measurements could constrain the positions of quasars to 
less than a milliarcsecond remains to be seen, but the VLBA 
typically measures the positions of sources about 10 microarc- 
seconds, so the SKA could in principle achieve 60 microseconds 
with its larger minimum wavelength and smaller size, detect- 
ing Earth-mass lenses out to about 10 pc with a source impact 
parameter of 0.1 milliarcseconds. If one were especially lucky 
and found an especially close encounter between the lens and 
source, the maximal displacement is 

in principle detectable with the SKA out to 250 pc. However, 
at such a distance the lens subtends such a small angle that 
finit e-sou r ce effec ts are likely to be important. 

iHevll (|2010bl l calculated the expected event rate of sub- 
stellar objects lensing bulge stars in the OGLE-II catalogue 
that can also be detected with the Square Kilometre Array 
(about 80,000 stars) under the assumption that the density of 
substellar objects in the disk of our Galaxy is about one-tenth 
of the total density. The total optical depth of such lenses 
is about 2 x 10~^. The calculation neglected the possibility 
of substellar objects in the Galactic halo and assumed that 
a lensing event lasts one day. Under these assumptions, the 
event rate where a source and lens align to within one Ein- 
stein radius is about once per 14 years. The lensing of bulge 
stars could yield an astrometric signature of diffractive mi- 
cr olensing in ad dition to the magnification signature discussed 
in lHevll(|2010bl ). 

This event rate can be scaled to the microlensing rate for 
fiat spectrum radio quasars for which the astrometric signature 
may be easier to measure. To achieve a rate of once per decade 
about 100,000 radio sources would have to be monitored. Be- 
cause one expects the lenses to be restricted to the plane of the 
Galax;y more or less, only those quasars that lie within ten de- 
grees of the Galactic equator should be considered as sources 
(about two steradians). The number counts of fiat-spectrum 
radio quasars giv e the expected event r ate as a function of the 
quasar radio fiux. lde Zotti et al] (;2010') compile the number of 
counts of radio-loud quasars at several frequencies. In partic- 
ular above several GHz the sample will be dominated by fiat- 



spectrum sources; furthermore, this is wher e the SKA is sen- 
sitive. From the number counts at 8.4 GHz jWindhorst et al.l 
ll993l : lFomalont et al. l l2002l : lHenkel fc Partridgell2005l ) one can 
estimate that there are about 10^ fiat-spectrum radio sources 
within ten degrees of the Galactic plane with fiuxes greater 
than 4 mjy. To get a sample of 10® sources requires a fiux 
limit of 0.3 mJy. What remains to be seen is how well future 
instruments will be able to centroid such sources as a function 
of fiux. These limits of course refer to the fiux of the entire 
source not just the high brightness temperature components 
that will show the most dramatic astrometric signatures. 



4 CONCLUSIONS 

The continuous monitoring of compact, distant radio sources 
may provide new way to probe the constituents of our so- 
lar neighbourhood, in particular freely, floating sub-stellar ob- 
jects. The astrometric signatures of diffractive microlensing 
can provide an estimate of the mass, distance and proper 
motion of the lensing object, possibly allowing follow-up ob- 
servations of the lens itself. Astrometric lensing even with- 
out diffraction effe cts can provide this information as well 
I Wambsgansd [2006f ) : however, diffraction typically amplifies 
the astrometric signature and radio observations often offer 
much higher angular resolution on the order of ten milliarc- 
seconds versus several hundred milliarcseconds in the optical. 

This letter has used the specifications of the SKA as a 
benchmark. Clearly the high angular resolution and high fre- 
quency offered by the SKA are helpful for the detection of as- 
trometric lensing in the radio; however, the high sensitivity of 
the SKA may not strictly be necessary if one focuses on bright 
radio sources. Perhaps, a purpose-built very-large baseline ar- 
ray of phased dipoles could achieve the needed angular resolu- 
tion (and possibly even a finer resolution than the SKA) with 
a sufficient sensitivity to continuously determine the centroids 
the brightest radio sources to the needed accuracy to detect 
low-mass objects in the solar neighbourhood. Furthermore, 
such a monitoring campaign could yield new insights on quasar 
physics as well as other ancillary results. The low expected op- 
tical depth for these events of about 2 x 10~^ would required 
the monitoring of 100,000 radio sources to achieve even the 
modest event rate of once per decade. These sources could be 
quasars or bulge giants, although the effect should be more 
pronounced with the high brightness-temperture quasars. 
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